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Sterile neutrinos of keV masses are one of the most promising candidates for the warm 
dark matter, which could solve the small-scale problems encountered in the scenario of 
cold dark matter. We present a detailed study of the production of such sterile neutrinos 
in a supernova core, and derive stringent bounds on the active-sterile neutrino mixing 
angles and sterile neutrino masses based on the standard energy-loss argument. 
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1. Motivation 

The observations of rotational curves of galaxies, the bullet cluster, gravitational 
leasing effects and cosmic microwave background have provided us with robust ev¬ 
idence that the matter content of our Universe is dominated by the non-baryonic 
dark matter. So far, a lot of attention has been focused on the scenario of cold dark 
matter (CDM), which has a negligible velocity dispersion at the radiation-matter 
equality and damps structures below the Earth-mass scales.— The favourite candi¬ 
dates for CDM stem from the well-motivated theories, which have been proposed to 
solve fundamental problems of the standard model of elementary particle physics,— 
such as the lightest supersymmetric particle and the axion. However, the CDM 
scenario suffers from a few serious problems in the galaxy and small-scale struc¬ 
ture formation, e.g., the overprediction of the observed satellites in the galaxy-scale 
halosi^ and the high concentration of dark matter in galaxies.— In the scenario of 
warm dark matter (WDM), a light-mass particle with a large velocity dispersion 
can suppress the structure formation up to the galaxy scales and thus solve the 
potential small-scale structure problems.— 

Sterile neutrinos of keV masses are a promising candidate for the WDM.— As 
a simple but instructive example, Dodelson and Widrow have proposed that right- 
handed neutrinos with masses ^ keV can be produced via neutrino oscillations 
in the early Universe and account for all the dark matter,— if they mix with the ordi¬ 
nary neutrinos via a tiny mixing angle 9 ^ in vacuum. Due to such a small 

mixing angle, sterile neutrinos can never be in thermal equilibrium. In the presence 
of a primordial lepton asymmetry, Shi and Fuller have observed that the production 
rate of sterile neutrinos could be enhanced by the Mikheyev-Smirnov-Wolfenstein 
(MSW) effect^i^ and the correct relic abundance can be obtained even for much 
smaller mixing angles.— One possible way to detect sterile-neutrino WDM is to 
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search for the X-rays from their radiative decays in the DM-dominated galaxies.— 
Conversely, the non-observation of an X-ray line from the local group dwarf galax¬ 
ies has placed restrictive limits on the mass and mixing angle of sterile neutrinos. 
Recently, the observation of an X-ray line around 3.5 keV in the Andromeda galaxy 
and the Perseus galaxy cluster has been claimed by two independent groups.— 
However, this claim is criticized by other authors, who attribute this X-ray line 
to the Potassium and Chlorine emission lines from the plasma in the galactic re¬ 
gions and the intergalactic gas.— At the present time, it is better to leave this 
problem open, and the future X-ray data could hopefully offer us a clue. Other 
limits can be obtained from the observations of Lyman-alpha forest and Supernova 
(SN) 1987A.— Put all together, the window for the Dodelson-Widrow mechanism of 
non-resonant production is closed, while the Shi-Fuller mechanism of resonant pro¬ 
duction is still viable.— Roughly speaking, sterile neutrinos with nis = 1 ^ 10 keV 
and 9 = could be WDM. However, it should be noticed that the observa¬ 

tional constraints depend crucially on the production mechanisms and evolution of 
sterile neutrinos in the early Universe, so the constraints can be evaded as in several 
new-physics models.— “— 

Moreover, the WDM sterile neutrinos could play an important role in generating 
the supernova asymmetries for pulsar kicks,— and in supporting supernova explo¬ 
sions.— >21 Hence it is interesting to reexamine the SN bounds on the keV-mass 
sterile neutrinos by studying in detail the production and propagation of sterile 
neutrinos in the SN core. On the other hand, as we show later, the phenomenon of 
neutrino oscillations and interactions in a dense medium is intriguing by itself. 

2. General Formalism 

The matter density in a SN core p = 3.0 x 10^^ g cm“^ is so high that even the 
weakly-interacting neutrinos could not escape freely.— Hence the production of 
sterile neutrinos in the dense matter can be very efficient through both neutrino 
oscillations and repeated scattering of ordinary neutrinos off background particles. 
In consideration of both coherent flavor oscillations and the decoherence caused by 
frequent scattering, it is convenient to describe the whole neutrino system in terms 
of density matrices.— We assume that neutrino interactions do not affect substan¬ 
tially the medium, whose configuration is mainly determined by the conditions of 
thermal equilibrium. The time duration of collisions between neutrinos and matter 
particles is much shorter than the evolution time scale, on which neutrino density 
matrices experience a significant change. On the other hand, the evolution of den¬ 
sity matrices is sufficiently rapid compared to the macroscopic and hydrodynamic 
time scale. 

Under the above assumptions, we follow the formalism mainly developed by 
Raffelt and Sigl,— and define the ensemble average of the density matrices 
of n-flavor neutrinos as (a](p)aj(q)) = (27r)^(5^(p - q)(pp)y and (5|(p)5j(q)) = 
(27r)3(53(p - q)(pp)y for i,j = 1,2, ••• ,n, where aj(p) and ^^(p) stand for the 
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annihilation operators for neutrinos and antineutrinos, respectively. The diagonal 
elements of and are just the occupation numbers, while the off-diagonal ele¬ 
ments encode the phase information. The evolution of the matrix of densities p^ is 
governed by— 


Pp = -i ^p] 


1 

+ ^ E / { [G>p^G“(l - pp) + h.c.] 

- [ppG“(l-pp,)G“ + h.c.] W“p,} . 


( 1 ) 


On the right-hand side, the first term contains neutrino kinetic energy and matter 
potential from the coherent forward scattering. The second term arises from the 
charged-current production and absorption of neutrinos, where is the projection 
matrix onto the neutrino flavor i, Vp the production rate, and Ap the absorption 
rate. The third term represents the contribution from scattering processes Pp -I- a —>■ 
i/p, +a' with a and a' being any kind of background particles, where G“ is a diagonal 
n X n matrix of couplings gf for the interaction type a and (pf)^Wp-p denotes the 
transition probability. In the case of antineutrinos, the evolution equation for Pp 
can be obtained in a similar way. 

Applying this general formalism to an ensemble of a keV-mass sterile neutrino 
mixing with only one flavor of ordinary neutrinos, we can simplify significantly the 
Boltzmann-type equations. To further simplify the picture, we consider a homo¬ 
geneous and isotropic SN core, in which neutrinos are essentially trapped and well 
in thermal equilibrium with the medium. Due to the isotropy of the medium and 
neutrino ensemble, the relevant interaction rates depend only on neutrino energies. 
The next important step is to average (p^;)^ and (p^:)^ over many cycles of flavor 
oscillations, since the oscillation length is much shorter than neutrino mean free 
path. More explicitly, we can estimate the oscillation length as 


Aosc ^ 0.7 cm 


E 


30 MeV / V sin26» 


10 


-4 


10 keV 


( 2 ) 


and the mean-free-path as 


-^mfp — 




1.1 X lO"^ cm 


30MeV\ 

Pl4 


E 


- 


( 3 ) 


where denotes the number of baryons, a 


uN 


GpE’^/ir the cross section of 


neutrino-nucleon interaction and Pi^ the matter density in units of 10^^ g cm“^. 
For a typical mixing angle sin^ 29 = 10“® and sterile neutrino mass nig = 10 keV, we 
always have ^ A^g^., corresponding to the weak-damping limit. This compar¬ 
ison indicates that neutrinos oscillate a large number of times before a subsequent 
scattering with nucleons. The averaged matrix of densities p^ can now be pa¬ 
rameterized in terms of the distribution functions of active and sterile neutrinos, 
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denoted as and /|;, respectively. If the active neutrino stays in thermal equilib¬ 
rium and if the mixing angle is so small that sterile neutrinos escape from the SN 
core immediately after production, the evolution equation of /J. is given by— 


fh = ^4,. 


E'^dE' 

27r^ 




( 4 ) 


where S 2 g = sin 20^ with 0^ being the neutrino mixing angle in matter, is 
the coupling constant for active neutrinos. Therefore, the emission rate of lepton 
number can be obtained by integrating the above equation over sterile neutrino 
energy. Similarly, the energy-loss rate can be calculated by convolving Eq. (4) with 
neutrino energy and integrating them over the whole energy range. 


3. Sterile neutrinos in SN cores 

Sterile neutrinos with masses in the keV range can be copiously produced in the SN 
core. For rris ^ 100 keV, the vacuum mixing angle of sterile neutrinos is stringently 
constrained sin^ 20 < 10“® in order to avoid excessive energy loss.— For smaller 
masses, however, the MSW effect on active-sterile neutrino mixing becomes very 
important and the SN bound on vacuum mixing angle is not that obvious. Note that 
the bounds on mixing angles depend on which neutrino species the sterile neutrino 
mixes with. In the following, we shall discuss and mixing cases, and 

emphasize the main difference between them. 


3.1. Ut-Vs mixing 

First, we concentrate on the SN bound in the simplest case of i^t-Es mixing, because 
Er and Er only have neutral-current interactions and essentially stay in thermal 
equilibrium with the ambient matter. In the weak-damping limit, which is always 
valid for supernova neutrinos mixing with keV-mass sterile neutrinos, the evolution 
of Et number density is determined by— 





E'‘^dE' 

2tt'^ 


We'eIe' 


( 5 ) 


where is the occupation number of Et^ and Y^e'e transition probability for 
e{E') + a ^ ^{E) + a with a being background particles in the SN core. Note that 
Eq. (5) is actually an immediate consequence of integrating Eq. (4) for Et-Es mixing 
over the sterile neutrino energy. For e^, the coupling constant is just g,- = 1, so we 
write the transition probability as Weie = yv^E'E- The interaction types include 
neutrino-nucleon and neutrino-electron scattering, whereas the former dominates 
over the latter. In a similar way, we can derive the evolution equation of the Ej- 
number density, involving the mixing angle the occupation number and the 
transition probability IF^/g. 
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Due to the MSW effect, the mixing angle of neutrinos in matter is different from 
that of antineutrinos, i.e.. 


sin^ 29u,v 


sin^ 26 

sin^ 26 + {cos2e ^ E/E,Y ’ 


( 6 ) 


where 6 denotes the vacuum mixing angle, and the sign “=p” refers to v and u. 
However, the actual sign depends on the lepton asymmetries in the matter, i.e., the 
matter potential The resonant energy = /S.m?/{2\Vu^\) can be written as 

E, = 3.25 MeV IFq - , (7) 

where Yq = (f - Ye - 2Y^,JjA, and 14 = {N^ - Ns:)/Nb with Nb being 
the baryon number density, and being the number densities of parti¬ 
cle X and its antiparticle x. As for tau neutrinos, the matter potential = 
— {Gf/V2)Nb {1 — Ye — 2Y^^ — 4Y^^) is negative if the typical values of Ye = 0.3, 
Y^^ = 0.07 and Y,y^ = 0 for a SN core are taken. Therefore, the mixing angle for 
Dr is enhanced by matter effects, and the emission rate for Dr exceeds that for i>r, 
indicating that a Vr-Dr asymmetry (i.e., Y^^ ^ 0) will be established. An interest¬ 
ing feedback effect emerges: (i) The chemical potential for tau neutrinos develops 
and thus changes the occupation numbers of Vr and Dr'., (ii) The Vr-Dr asymmetry 
shifts the resonant energy and thus modifies the mixing angles 6i, and (iii) 
Both effects in (i) and (ii) will feed back on the emission rates. Hence a stationary 
state of this active-sterile neutrino system could be achieved if the emission rates for 
neutrinos and antineutrinos become equal to each other.— However, whether such 
a stationary state can be really reached depends on the mixing angle and sterile 
neutrino mass. In any case, it is straightforward to calculate the energy-loss rate 
by following the time evolution of the neutrino system. 


3.2. v'e-v's mixing 

Since there is a great similarity between muon and tauon neutrino interactions in a 
SN core, we expect that mixing case shares the same features as the previous 
case. Then, we turn to the Ue-Vg mixing case, where the physical processes for the 
production and evolution of sterile neutrinos are quite different. 

First, the charged-current interaction e~ +p ^ Ve+n leads to a direct production 
of r'e, and an absorption of Ue by the medium as well. In addition to the charged- 
current interaction, neutrino scatterings off degenerate electrons and non-degenerate 
nucleons are also important. Therefore, both terms on the right-hand side of Eq. (4) 
exist for the case of Ve-i's mixing. 

Second, the emission of Vg reduces the electron lepton number in the SN core, 
which is in thermal and beta equilibrium. From the latter condition, we can get a 
relationship pe = + A among the electron chemical potential pe, the electron- 

neutrino chemical potential and the difference A = 6n — Dp between the neutron 
fin and proton fip chemical potentials. The parameter A = 50 ^ 100 MeV is 
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determinable from the equation of state of the nuclear matter. In our discussions, 
we assume the beta equilibrium is always guaranteed and the equation of state is not 
significantly modified by the emission of sterile neutrinos. Under these assumptions, 
/t = 55 MeV is taken to be valid throughout. The emission of sterile neutrinos 
affects the number density of and thus and which in turn modify the 
matter potential for neutrinos and the mixing angle in matter. The modification of 
neutrino mixing angle obviously feeds back on the emission rate. 

Third, since and i>t are always generated in pairs, there is no asymmetry 
between them, namely = 0, and likewise for and i>^. As a result, the 
resonant energy = Am^/(2||) for neutrino oscillations is found to be 



where = 0 is implemented. For the initial conditions /ie = 218 MeV 

and = 163 MeV, one obtains Y^ = 0.3 and Y^,^ = 0.07 for the temperature 
T = 30 MeV and matter density p = 3.0 x 10^^ g cm“^. Unlike in the Vt-Vs mixing 
case, the matter potential 1/^^ is positive, implying the resonance is initially in the 
neutrino sector. Therefore, the emission rate of Vs is at the beginning much larger 
than that of Vg, and the lepton number decreases rapidly. But even small changes 
of Yg and are able to drive to infinity, where 1 — 3Yg — 4Yj^^ = 0 is reached. 

It is worthwhile to mention that a large chemical potential indicates the 
population of bg is rare. The initial emission of Ug and the corresponding energy 
loss in this period become crucial for us to draw any SN limits based on the energy- 
loss argument. When is approaching infinity, the mixing angles in matter for 
neutrinos and antineutrinos are equal, and both of them are close to the vacuum 
mixing angle. During the evolution, the matter potential flips its sign and the 
resonance moves to the antineutrino sector. However, the bg number density is 
largely suppressed due to the large chemical potential, even the enhancement of 
antineutrino mixing angle cannot enlarge the total emission rate remarkably. 

4. SN bound on sterile neutrinos 

Given the sterile neutrino mass mg and vacuum mixing angle 0, the energy loss rate 
£{t) due to sterile neutrino emission can be calculated by following the evolution 
of Vt-Vt asymmetry Yi,^(t). It has been found that the stationary state can be 
reached within one second and the feedback effect is very important for 20 keV < 
mg < 80 keV and 10“® < sin^ 29 < 10“^. To avoid excessive energy losses, we 
require that the average energy-loss rate {£) = /J"'* £ (t) dt with Td = I s should be 
{£) < 3.0 X 10^^ erg cm“^ s“^. Otherwise, the duration of neutrino burst from SN 
1987A would have been significantly reduced. 

In Fig. 1, we show the contours of energy-loss rates in the (sin^ 20, mg)-plane, 
where we have assumed a homogeneous and isotropic core with matter density 
p = 3.0 X 10^'^ g cm“^ and temperature T = 30 MeV. Based on the energy-loss 
argument, the shaded purple region has been excluded. The most stringent bound 
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IQ-IO jq-8 jq-6 jq-4 jq-2 1 

sin^ (2 6 ) 


Fig. 1. Supernova bound on sterile neutrino masses rris and active-sterile mixing angles 9 in the 
Ur-^s mixing case, where the shaded purple region is excluded by the energy-loss argument while 
the shaded green one by the energy-transfer argument™. The excluded region will be extended 
to the dashed (red) line if the build-up of degeneracy parameter is ignored, i.e., 77 (f) = 0. The 
dot-dashed (green) line represents the keV-mass sterile neutrinos as dark matter with the correct 
relic abundance = 0 . 1 . 

sin^ 26 < 10“® arises for nis = 50 keV. For the large-mixing angle region, the 
energy-loss rate is actually small, because sterile neutrinos have been trapped in 
the core and cannot carry energies away. However, the mean free path of sterile 
neutrinos is comparable to or even larger than that of ordinary neutrinos, indicating 
that they may transfer energies in a more efficient way. As a consequence, the 
duration of neutrino burst will be shortened by emitting neutrinos more rapidly. 
In this sense, the excessive energy transfer should be as dangerous as the excessive 
energy loss. Therefore, the large-mixing angle region is excluded when the energy- 
transfer argument is applied. The green dot-dashed line in Fig. 1 indicates the 
relic abundance of dark matter = 0.1, where keV-mass sterile neutrinos are 
warm dark matter and the non-resonant production mechanism is assumed. If we 
ignore the feedback effect (i.e., a vanishing chemical potential for tau neutrinos 77 = 
= 0 ), the excluded region will extend to the red dashed line, which overlaps 
the relic-abundance line. However, the mixing angles are essentially unconstrained 
in the favored warm-dark-matter mass range 1 keV ^ mg < 10 keV. 

As for the j^g-j^s-mixing case, we have included all the scattering processes of 
neutrinos, and solved the Boltzmann-like equation of the occupation numbers. In 
assumption of thermal equilibrium of Vg, one can derive the evolution equation 
of the degeneracy parameter 77 = It is straightforward to verify that the 

infinite is reached quickly, and the crossing from a negative to a positive matter 
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> 10 

1 


10 '^^ 10''“ 10'^ 10'^ lO'"^ 10'^ 

sin^ (29) 



Fig. 2. Supernova bound on sterile neutrino masses rus and active-sterile mixing angles 8 in the 
mixing case. Based on the energy-loss argument, the parameter space between two dot- 
dashed curves is excluded. The region of large mixing angles at the upper-right corner is also 
disfavored if the energy-transfer argument is implemented. 


potential takes place. In Fig. 2, using the same method as before, we have presented 
the contour plot of the energy-loss rate in the parameter plane of sterile neutrino 
masses and active-sterile neutrino mixing angles.— On the upper-right corner, the 
energy-transfer argument can be applied to exclude this large-mixing angle region. 
For sterile neutrinos of masses above 1 keV, the leftmost exclusion line is nearly 
vertical around sin^ 26 = 10“®. For much smaller masses mg <C 1 keV, the emission 
rate is highly suppressed by the matter effect, so the bound is rather loose. 

5. Summary and Outlook 

Since sterile neutrinos of keV masses as a promising candidate for warm dark matter 
are interesting in astrophysics and cosmology, we have considered their production 
in a SN core, and derived restrictive bounds on their masses and mixing angles by 
requiring the excessive energy loss caused by sterile neutrinos to be under control. 

The main idea is to follow the time evolution of occupation numbers of neutrinos, 
and compute the emission rates of lepton number and sterile neutrino energy. Our 
final results about SN bounds are presented in Fig. 1 and Fig. 2. We have also 
emphasized the feedback effects in both Vt-Vs and mixing cases: The emission 
of sterile neutrinos leads to changes in the lepton number asymmetries and thus the 
neutrino matter potentials, which in turn affect the emission rates via the MSW 
effects on neutrino mixing angles. In the case of mixing, the MSW resonance 
occurs in the antineutrino sector, and thus an asymmetry between and Dr number 
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densities will be developed. This asymmetry contributes to the matter potential, 
intending to modify the mixing angles in matter and equating the emission rates of 
Vs and Vs- In the case of v^-Vs mixing, the matter potential is rapidly driven to 
zero, implying that the neutrino and antineutrino mixing angles are close to the one 
in vacuum, which is in agreement with the result found in the literature.— However, 
for much smaller masses, the emission rates will be suppressed by the matter effects, 
and the SN bound becomes very weak. 

It is worth pointing out that keV-mass sterile neutrinos could also be produced 
in the collapsing phase of massive stars.— Such sterile neutrinos after production 
propagate from the inner core to the outer one, and oscillate into electron neutrinos, 
which could deposit their energies in the out-layer matter and support SN explo¬ 
sions. A dedicated study of keV-mass sterile neutrinos in both core-collapse and 
cooling phases is certainly necessary and intriguing, and left for future works. 
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